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Abstract 

We systematically investigate the topological band structures of half-Heusler compounds using 
first-principles calculations. The modified Becke-Johnson exchange potential together with local 
density approximation for the correlation potential (MBJLDA) has been used here to obtain accu- 
rate band inversion strength and band order. Our results show that a large number of half-Heusler 
compounds are candidates for three-dimensional topological insulators. The difference between 
band structures obtained using the local density approximation (LDA) and MBJLDA potential is 
also discussed. 



PACS numbers: 71.15.Mb, 71.20.Nr, 73.20.At 



I. INTRODUCTION 



Motivated by their potential applications in spintronics and quantum computing, the 
search for three-dimensional topological insulators (3DTI) has attracted considerable the- 
oretical and experimental intereslP^. These materials are called "topological" because 
they are distinguished from ordinary insulators by the so-called Z 2 topological invariants 
associated with the bulk band structur e 1 14 * 15 ! On the theory side, the state-of-the-art first- 
principles calculations guided by topological band theory^HSl provided a powerful tool 
for uncovering new families of 3DTI. Based on this approach, it has been recently predicted 
that ternary half-Heusler compounds can realize the topological insulating phase with proper 
strain engineering^^. 

There are several ways to determine the band topology of an insulator. Intuitively, one 
can count the number of band inversions within the entire Brillouin zone - an odd number 
indicates that the material may be a 3DTPS1. This method depends on an accurate interpre- 
tation of the atomic origin of the bands and is better suited for crystals with a high-symmetry 
lattice. A more rigorous method is to directly evaluate the Z 2 topological invariants. For 
materials with inversion symmetry, the parity criteria developed by Fu and Kane^ can be 
readily applied. On the other hand, if the inversion symmetry is absent, one must resort to 
the more elaborated lattice computation of the Z 2 invariants^. We have recently applied 
this method to study the distorted half-Heuslers^ and the noncentrosymmetric chalcopy- 
rite compounds^. In practical calculations, both methods require an accurate knowledge of 
the bulk band structures. So far, previous works on 3DTI have employed either the local 
density approximation (LDA)PH22] or generalized gradient approximation (GGA)P^ for the 
exchange-correlation potential. However, it is well known that these approximations have the 
tendency to underestimate the bandgap. In particular, if the material has a small positive 
bandgap, this underestimation may yield a negative value. Therefore, these approximations 
may falsely predict an inverted band structure when the band order is actually normal. 

Recently, a new semilocal potential that combines the modified Becke- Johnson (MBJ) 
exchange potential and the LDA correlation potential, called MBJLDA^, was proposed to 
obtain accurate bandgaps and band order. The MBJLDA potential are computationally as 
cheap as LDA or GGA, but it has similar precision compared with the more expensive hybrid 
functionals and GW method. It has been demonstrated 2 - 4 that the MBJLDA potential can be 
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used to describe many types of solids, including wide bandgap insulators, small bandgap sp 
semiconductors, and strongly correlated 3d transition-metal oxides. More importantly, the 
MB JLDA potential can effectively mimic the behavior of orbital-dependent potential around 
the bandgap, so it is expected to obtain accurate positions of states near the band-edge, 
which are the keys to determine the band inversion and the band topology 

In this work we perform a systematic investigation of the band topology of the half-Heusler 
family using the MB JLDA potential, and compare it with the LDA result. The improved 
accuracy of the MBJLDA potential over LDA is tested in the ordinary insulator CdTe and 
the topologically nontrivial 3D-HgTe, whose band structures are known experimentally and 
can be used as a benchmark. We then focus on the ternary half-Heusler compounds. By 
calculating the band inversion strength, we confirm our previous prediction that a large 
number of half-Heuslers are possible candidates for 3DTI or 3D topological metaP^. The 
difference between LDA and the MBJLDA potential is also discussed and a clear discrepancy 
is found between these two methods. The sensitivity of the calculated band structure due 
to the type of exchange-correlation potential calls for more detailed experimental works. 

II. METHOD 

The band structure calculations in this work were performed using full-potential linearized 
augmented plane-wave (FP-LAPW) method^, implemented in the package wien2kP^. A 
converged ground state was obtained using 10,000 k points in the first Brillouin zone with 
K max RMT = 9.0, where Rmt represents the muffin-tin radius and K max is the maximum 
size of reciprocal-lattice vectors. Wave functions and potentials inside the atomic sphere are 
expanded in spherical harmonics up to / = 10 and 4, respectively. Spin-orbit coupling is 
included by a second-variational procedure^, where states up to 9 Ry above fermi level are 
included in the basis expansion, and the relativistic p±/2 corrections were also considered for 
hp and 6p orbit in order to improve the accuracj^Zl. 
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(a) CdTe-LDA (b) (c) CdTe - MBJLDA 




FIG. 1: (color online)Band structure of CdTe and HgTe. (a) LDA and (c) MBJLDA results of 
CdTe; (d) LDA and (f) MBJLDA results of HgTe. The r 6 , T 7 , and T 8 state are denoted by 
red, black, and blue color, respectively. The size of dots is proportional to the probability of s- 
orbit projection, (b) and (e) is schematically experimental band structure of CdTe and HgTe, the 
red solid dots denote the s-orbit originated Tq state. The band inversion strength is defined as 
AE = Er 6 — Er 8 - The calculated and experimental energy of Fq, Fj, and r 8 state are listed in 
Table ffl 

III. RESULTS AND DISCUSSIONS 
A. CdTe and HgTe 

We first test the MBJLDA potential for the binary compounds CdTe and HgTe. The 
energy bands of both compounds at the T point near the Fermi level split into T 6 (2-fold 
degenerate), (2-fold degenerate), and r 8 (4-fold degenerate) states due to the zinc-blende 
crystal symmetry and strong spin-orbit interaction. The experimental band order of CdTe 
and HgTe are (from high to low energy) T 6 , r 8 , T 7 (Ref. 1281291) and r 8 , r 6 , T 7 (Ref. I30j) . 
respectively. They are schematically shown in Fig. [ijb) and[lje). From the viewpoint of 
band topology, the main difference between these two compounds is that CdTe possesses 
a normal band order, i.e., the s-like r 6 state sits above the p-like T 8 state, while HgTe 
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TABLE I: The energy of T 6 , T 7 , and T 8 state of CdTe and HgTe calculated with LDA and MBJLDA 



at their experimental lattice constants 6.48Aand 6.46A, respectively^. Experimental values are 
also shown for comparison. The fermi level are always located on T% state, and £"r 8 are set to zero. 
The energy units are eV. 
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possesses an inverted band order in which the T 6 state is occupied and sits below the T 8 
state. We then define the band inversion strength AE as the energy differences between 
these two states, i.e. 

AE = E F6 - E Ts . (1) 

A negative AE typically indicates that the materials are in a topologically nontrival phase, 
while those with a positive AE are in a topologically trivial phase. 

Figure [T] shows the LDA and MBJLDA band structures of CdTe and HgTe at their 
experimental lattice constants 6.48Aand 6.46A, respectively^. The calculated energy of the 
Tq, Tj, and Tg states together with their experimental values are listed in Table [Tj For CdTe, 
the LDA potential yields a serious underestimation of AE. For HgTe, it obtains the wrong 
band order of the and states (but it does predict an inverted band order). On the 
other hand, the MBJLDA result shows an excellent agreement with experiments for both 
compounds. We therefore conclude that the MBDLDA potential is better suited to calculate 
the topological band structure. 
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FIG. 2: (color online) Crystal structure of half-Heusler compound XYZ in the F43m space group. 
Green spheres at (0.5,0.5,0.5) are atom X, dark blue spheres at (0.25,0.25,0.25) are atom Y, and 
pink spheres at (0,0,0) are atom Z. 

B. Half-Heuslers 

Having established the improved accuracy of the MBJLDA potential over LDA, we now 
turn to the ternary half-Heusler compounds described by space group FA3m. The chemical 
formula of these materials is XYZ, where X and Y are transition or rare earth metals 
and Z a heavy element. It can be regarded as a hybrid compound of XZ with rock- 
salt structure, and XY and YZ with the zinc-blende structure (see Fig. [2|. In our band 
structure calculations, the lattice constants are taken from experimental data library^ESl. 
For compounds without experimental lattice constant, we use the value obtained by total 
energy minimization in first-principles calculations. 

The band structures of half-Heuslers are very similar to CdTe/HgTe. In particular, the 
low-energy electron dynamics is dominated by energy bands at the T point. Therefore, 
the band topology of half-Heuslers can be characterized by the band inversion strength 
AE defined in Eq. ([!]). In order to systemically explore their topological phase, we have 
calculated AE by both LDA and MBJLDA for 24 half-Heusler compounds, as shown in Fig. 
[3]^a) and (b), respectively. Following our previous work (Ref. ITT]) , we have also calculated the 
Z 2 invariant for all half-Heusler compounds investigated here. The calculated Z 2 invariant 
agrees with the intuitive band inversion picture, that is, materials with negative AE are 
in topologically nontrivial phase. As shown in Fig. [3|b), we have identified 9 half-Heusler 
compounds, LuPtSb, LaPtBi, YPtBi, ScPtBi, ThPtPb, LaAuPb, YAuPb, LuPdBi, LuPtBi 
as possible candidates for topological insulators or topological metal. (Note that ScPtBi, 
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FIG. 3: Band inversion strength of half-Heusler compounds calculated by LDA (a) and MBJLDA 
(b). The materials with open symbols indicate that there is not experimental reports, and the 
lattice constant are obtained by first-principles total energy minimization. The topological insulator 
or metal (after applying a proper uniaxial strain) with AE < are located below the horizontal 
line; and the ordinary insulator or metal with AE > are located above the horizontal line. 
The squares mark insulator or zero-bandgap semiconductor and diamonds mark the metal. The 
materials enclosed by circle indicate that their topological phase predicted change by using different 
exchange-correlation potential. 

ThPtPb, LaAuPb here are virtual compounds.) Here, the band structures of the first 6 
compounds are with zero band-gap similar to that of HgTe, and the last 3 compounds 
possess the nontrivial metallic band structures similar to Fig. He) or the conduction bands 
droping below the Fermi level near the X point. 

As shown in section III. A, the LDA result generally yields a smaller value of the band 
inversion strength AE when compared with the MBJLDA result, which agrees much better 
with the experimental data. We emphasize that, unlike the usual energy gaps, here AE 
can take on both positive and negative values, hence a positive AE is always larger than 
a negative one. This trend is also confirmed for half-Heuslers. As shown in Fig. [3j AE by 
MBJLDA is always larger than that by LDA. In particular, if LDA calculation predicts a 
small negative AE, it may become positive when using MBJLDA potential, such as ScAuPb, 
and YPdBi. This suggests that when predicting topological phases, one must be careful with 
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(e) YPtBi-LDA (f) YPtBi - MBJLDA (g) LaPdBi - LDA (h) LaPdBi - MBJLDA 




FIG. 4: (color online) LDA and MBJLDA band structure of half-Heusler compounds LaPtBi(a)(b), 
YPdBi(c)(d), YPtBi(e)(f), and LaPdBi(g)(h). The labels are the same as Fig. [lj 

the types of exchange-correlation potentials. 

To illustrate the general difference between the LDA and MBJLDA results, we present 
four scenarios shown in Fig. [4] using some typical examples. Both LDA and MBJLDA 
calculations predict that LaPtBi [Fig. ^ a) and |4^b)] has an inverted band structure, which 
is very similar to HgTe. The main difference is the band order exchange between Tq and 
state. Although the band topology will not change, the LDA result clearly underestimates 
the AE due to the downward movement of r 6 state. For YPdBi, the problem is much 
more severe as the two potentials give different band topology. As shown in Fig. |4|c) and 
|4^d), the LDA result shows an inverted band order but the MBJLDA result shows a normal 
band order. The two types of exchange-correlation potentials also affects the electronic 
behavior around the Fermi level as shown for YPtBi in Fig. |4](e) and|4]^f). The metallic 
phase predicted by LDA calculation becomes a zero-bandgap semiconductor in MBJLDA 
calculation. The band order of T 6 and T 7 is also exchanged, just like LaPtBi. 
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FIG. 5: (color online) Band structures of LaPdBi by MBJLDA (a) without and (b) with a 7% 
hydrostatic strain. The application of a hydrostatic strain causes the s-like Tq bands (marked by 
red dots) to jump below the valence band top, providing the necessary band inversion that leads 
to a nontrivial topological order. The other labels are the same as Fig. [T] 

As mentioned in our previous workP^, for half-Heusler compounds with small band gaps, 
the nontrivial topologically phase can be generally realized by applying hydrostatic strain 
to change the band order. For example, LaPdBi has normal band order in its native state, 
a 7% change in the lattice constant converts the trivial topological phase into a nontrivial 
topological phase (Fig. [5]). Again, there are different effects on the band order under the same 
hydrostatic strain when using different exchange-correlation potentials. For example, when 
one stretches lattice constant to a = + 3%ao, LDA calculation will yield an inverted band 
structure [Fig. Qg)], while it still has a normal band order by MBJLDA calculation [Fig. 

m 



IV. SUMMARY 



In summary, we have systematically investigated the topological band structure of the 
half-Heusler family using both LDA and MBJLDA exchange-correlations potential. Our 
result shows that a large number of half-Heusler compounds are candidates for three- 
dimensional topological insulators. We also discussed the main differences between the 
LDA and MBJLDA potentials 
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